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While standard heavy fermion metals feature a single spin-1/2 local moment per unit cell, more 
complicated systems with multiple distinct local moments have been synthesized as well, with 
Ce3Pd2o(Si,Ge)6 being one example. Here, we discuss the physics of a Kondo lattice model with 
two local-moment sublattices, coupled with different Kondo couplings to conduction electrons. The 
phase diagram will be strongly modified from that of the standard Kondo lattice if the characteristic 
screening temperatures of the distinct moments are well separated. Therefore, we investigate the 
interplay between the two Kondo effects using a local self-energy approximation via slave bosons. 
We find that the two Kondo effects can either compete or co-operate depending on the conduction- 
band filling. In the regime of competition, small differences in the two Kondo couplings can lead to 
huge differences in the respective Kondo scales, due to non-trivial many-body effects. We also study 
the low-temperature properties of the collective heavy Fermi-liquid state and propose a connection 
to depleted Kondo lattice systems. 

PACS numbers: 71.27.+a, 72.15.Qm, 75.20.Hr, 75.30.Mb 



I. INTRODUCTION 

Inter-metallic compounds containing rare-earth ions 
have been the subject of intense experimental and the- 
oretical efforts, owing to their fascinating physical prop- 
erties. These include heavy fermion mass enhancement, 
unconventional magnetism and superconductivity, non- 
Fermi liquid behavior as well as the exciting possibility of 
quantum criticality beyond the Landau-Ginzburg- Wilson 
paradigm.EJEl Microscopically, most systems contain one 
rare-earth ion per unit cell, with a single unpaired / elec- 
tron or hole, lending themselves to a theoretical descrip- 
tion in terms of a Kondo lattice model with one spin-1/2 
local moment per unit cell. 

Recent worU^"SI on the family of ternary compounds 
(RE)3Pd2oX6 (Re=rare-earth, X=Si, Ge), which feature 
two inequivalent Kondo sites per unit cell,13i5i uncovered 
a variety of interesting phenomena. Upon cooling, two 
magnetic phase transitions occur .fiSHn] Ce3Pd2oSi6, 
regarded as one of the heaviest electron Kondo sys- 
tems with a low-temperature specific-heat coefficient of 
C/T = 8 J/mol Ce K^, the lower transition can be tuned 
to zero temperature by application of a magnetic field, re- 
sulting in non- Fermi liquid signatures with a linear tem- 
perature dependen ce of the resistivity.E^ISl Furthermore, 
substitution studied^'i^ suggest that quantum criticality 
might also be achieved by pressure. Together, these re- 
sults call for theoretical investigations of Kondo lattices 
with multiple Kondo ions. 

A. General considerations 

Theoretically, little is known about the behavior of 
Kondo lattices with multiple distinct local moments per 
unit cell. Here we start with a few general considerations 



for a system with two distinct Kondo ions. 

For sufficiently strong Kondo couplings, a paramag- 
netic Fermi liquid phase with fully screened moments 
can be expected. Its Fermi volume is given by Vfl = 
Kdintot mod 2) with ntot = 2nc + = 2nc + 2, where 
2nc and 2nf denote the average number of c and / 
electrons per unit cell, respectively. Further, Kd — 
{2tt)'^ / {2vo) is a phase space factor, with vq the unit cell 
volume and the factor of two accounting for the spin de- 
generacy of the bands. If, in contrast, inter-moment in- 
teractions dominate, phases with magnetic long-range or- 
der (LRO, e.g. antiferromagnetic) without Kondo screen- 
ing appear natural; alternatively, inter-moment singlet 
formation of spin-Pcierls type can also occur. Given two 
spin-1/2 moments per crystallographic unit cell (more 
generally, an even number), both phases can be realized 
without further lattice symmetry breaking, resulting in 
a Fermi volume Vaf — Kd{2ncinod2) which equals Vfl- 

Interestingly, a novel regime with partial Kondo screen- 
ing may be reached when the characteristic Kondo scales 
of the two ions are distinctly different. First, upon in- 
creasing the temperature in the paramagnetic Fermi liq- 
uid, the magnetic entropy will be released in two stages, 
allowing to define two "screening" temperatures Ti and 
T2 associated with the two ions. Second, upon increasing 
inter-moment interactions at low temperature, magnet- 
ically ordered phases can occur, where one set of mo- 
ments is strongly polarized with little screening, while 
the second set of moments still forms heavy quasipar- 
ticles with weak polarization. Such a magnetic phase 
is - apart from the possibility of Lifshitz transitions as- 
sociated with changes in the Fermi-surface topology - 
continuously connected to a local-moment metallic mag- 
net without Kondo screening, but observables will re- 
flect the co-existence of local-moment magnetism with 
heavy fermion behavior. A schematic phase diagram 



2 















^ ^ ^ 
partial screening 


/\ 


^coh / 


\ 

\ "light" LRO 


FL ~~~ ( f 


1 

■■ ► 



"heavy" LRO 

FIG. 1: Schematic phase diagram of a Kondo lattice with 
two inequivalent local moments per unit cell, as function of 
temperature T and the magnitude of the inter-moment ex- 
change I (measured relative to a microscopic Kondo scale, 
e.g., the largest single-ion Kondo temperature of the system). 
The crossover temperatures T\,2 mark the high-temperature 
onset of screening of the two ions in their lattice environment, 
while Tcoh is the lattice coherence scale below which the sys- 
tem displays Fermi-liquid behavior. Partial Kondo screening 
is realized for T2 < T < Ti (shaded). Upon increasing /, the 
system enters a long-range-ordered (e.g. antiferromagnetic) 
phase where the order is mainly carried by the Iow-Tk ions 
(here "2"), while the high-TK ions ("1") remain almost fully 
screened - this is a "heavy LRO" phase. Further increasing 
/ induces a crossover into a more conventional regime where 
all moments participate in the order ("light LRO"). It is con- 
ceivable that the "light" and "heavy" LRO regimes are char- 
acterized by different ordering patterns, necessitating further 
phase transitions inside the LRO phase. The focus of this 
paper is on the regime of small /. 



resulting from this discussion, which extends Doniach's 
consideration^^ to the case of multiple distinct Kondo 
ions, is shown in Fig. [T] 

A central question is then under which conditions the 
two screening scales Ti and T2 will be vastly different. 
On the one hand, this of course depends on bare model 
parameters (i.e. the Kondo couplings and local densities 
of states), but - on the other hand - there is room for mu- 
tual influence and non-trivial interplay of the screening 
processes. The latter can be phrased into the following 
question: Does the onset of Kondo screening for one set 
of moments at some temperature Ti promote or inhibit 
the screening of the second set of moments? Naive argu- 
ments can be found for either scenario: If the main effect 
of the higher-temperature Kondo screening is to inject 
additional states near the Fermi level, then the lower- 
temperature Kondo effect should be boosted. Alterna- 
tively, a variant of the exhaustion argumentfi^l would im- 
ply that, once the higher-temperature Kondo effect sets 
in, fewer conduction electrons are left available for the 
lower-temperature Kondo effect, which is then inhibited. 
In the latter scenario, small differences in the two Kondo 



couplings can cause largely different Ti_2. Naturally, the 
lower of the two screening temperatures Ti_2 sets an up- 
per bound to the low-temperature coherence scale Tcoh 
of the heavy Fermi liquid (experimentally often defined 
as the upper boundary of behavior in the resistivity) . 
However, even for standard Kondo lattices with one lo- 
cal moment per unit cell it has been showrP2H22l that T^oh 
can be tiny compared to the onset temperature of screen- 



ing (provided that inter-site effects are negligible^, in 
particular in the so called "exhaustion regime" of small 
band filling. Clearly, the issue of what determines the 
coherence scale Tcoh m the case of multiple distinct local 
moments is far from obvious. 

The purpose of this paper is to tackle the simplest 
set of questions outlined above, namely that of screening 
and coherence scales in the heavy Fermi-liquid regime 
of a Kondo lattice model with two distinct local mo- 
ments. These issues can be efficiently addresse d in a lo- 
cal self-energy approximation using slave bosons.l^^ We 
shall discuss the interplay between the Kondo crossovers 
of the two sets of moments. As detailed below, we 
find that - depending on the conduction-band filling 
ric - the Kondo processes may either compete or co- 
operate. For both cases, we describe the consequences 
for the low-temperature behavior and discuss experimen- 
tal implications. In particular, we shall argue that small 
conduction-band filling tends to open a window of par- 
tial Kondo screening as in Fig. [T] Further, the dominant 
influence of the conduction-band filling furnishes a link 
between the problem studied here and a seemingly differ- 
ent one, namely a Kondo lattice with random depletion 
of local moments.'^StSll 

We note that partial Kondo screening, somewhat simi- 
lar to the regime advertised in Fig. [T] has been predicted 
to occur in certain frustrated Kondo lattices, where a 
spontaneous modulation of the Kondo effect is instru- 
mental in releasing geometric frustration.l2Sl29| 



B. Outline 



The remainder of this paper is organized as follows. In 
Sec. |ll]we introduce the microscopic Kondo lattice model 
with two distinct Kondo ions, to be studied in this pa- 
per. This model is studied in the framework of the slave- 
boson mean-fleld approach. The emergent Kondo scales 
are derived and discussed in Sec. |III[ including the phys- 
ical picture of competing vs. co-operating Kondo effects. 
The low-temperature properties of the Fermi liquid are 
the subject of Sec. 



rV| We shall conclude by relating 
our results to depleted-Kondo-lattice physics and by dis- 
cussing directions for future work. Technical details are 
relegated to the appendices. 
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FIG. 2: Sketch of a lattice described by the model ([T]). Two 
types of local moments are coupled to the conduction elec- 
trons of the two sublattices with Kondo couplings Ji,2- The 
shaded area defines the unit cell including two spin- 1/2 local 
moments. Dots symbolize translation symmetry. 



II. MODEL AND MEAN-FIELD 
APPROXIMATION 

In the body of this paper, we discuss a Kondo lattice 
model with two inequivalent local-moment sublattices, 
with a Hamiltonian of the form 

% — — ^ ^ijcJ(TCjT+^ JiSi- Si — ^"^2 (^cj^Cia — ric^ , 

(1) 

in standard notation, with = X^o-o-' cJo-^o-cr'Cicr'/^ being 
the spin density of conduction electrons at site i and ric 
the average number of c electrons per site. The lattice, 
with a total number of unit cells M , is assumed to be 
bi-partite with sublattices A and B. The Kondo coupling 
Ji differs for both sublattices and is given by 

Ji = Jl^iA + J2kB (2) 

where 8ia — 1 (0) for i G a {i ^ a) and a — A,B. In the 
following, we will assume Ji > J2. 

The unit cell of the model ([T]) contains two c orbitals 
and two local-moment (or /) orbitals. Fig. [2j For most 
of the paper, we shall consider simple nearest-neighbor 
hopping only, such that = t on nearest-neighbor bonds 
and tij = otherwise. Both longer-range hopping and 
additional kinetic-energy modulations (compatible with 
the lattice symmetry) are not expected to qualitatively 
change our conclusions; this is explicitly shown in Ap- 
pendix |B] for the case of a square lattice with second- 
neighbor hopping. For Ji = J2 our model reduces to the 
standard Kondo lattice with conduction electron filling 
ric, dubbed "homogeneous" case in the following. 



A. Slave-boson mean-field approximation 

The model ([T]) will be tackled using a standard slave- 
boson technique. We use a fermionic representation for 
the local spin Si , introducing auxiliary fermions fia such 



that 

— ^ ^ /icr'^'o'CT' /ict' /2j 

together with a local constraint fiafia = 1. The 
Kondo interaction in ([T]) now takes the form 

= - ^ ^ J^flc^acl^,f,a', (3) 

where additional terms have can be absorbed into chem- 
ical potential and have been thus dropped. The interac- 
tion ^ is then decoupled using a field bi conjugate to 
(Ecr V'Jifjcr'^icr)- ^hc Lagrauglan then reads 

/: = T ^ (-(iw + fl)S,j ~ Uj) c\^Cj„ - ^ \^ftf^a 
iJ.LJ,a i.a 

(4) 

Recall that TV is the number of two-site unit cells of the 
lattice. 

The mean-field theory is obtained from the Lagrangian 
Q by assuming that the fields b and A have uniform and 
static expectation values on each sublattice 

(bia) = ba, (Xia) — Aq. 

The approximation ([5| is formally exact in the N ^ oo 
limit of a model with the spin symmetry extended to 
SU(iV). 

Performing a Fourier transformation, Cko-Q(/ko-a) — 
^ Ljsa e*'' '''9o-(/j<T), k running over the Brillouin 
zone of the Bravais lattice, one obtains the following 
mean-field Lagrangian 

^—lE E (i./M«.-i(»,k)(t- ) 

a,^ k,q,LJ,(7 

+ 2A/'^62 ^AA^A„ + 2A/'ne/x, (5) 

a a 

where the 4x4 inverse Green's function matrix is 
defined by 

^-1/- ^ f (i^ + Sap - SapCk) -y/Jaba Sal3 \ 

^"^^ ' ^ ^ -VJ^baSap {toj + Xa) Safi J ' 

(6) 

and 

is the kinetic-energy matrix whose diagonal (off-diagonal) 
elements result from intra- (inter-) sublattice hopping. 
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The mean-field free energy is readily obtained from ([s]) 
-rTrlog(-g-i)+2AA^&2 



MF 



(8) 



where Tr is a trace over all internal indices a, k, w, a. 

Self-consistent mean-field equations are obtained by 
minimizing the free energy ([s]) with respect to Xi, 
and /X. The result is: 



1 = 



2T 



k,a; 



(9) 



where the average occupation number of the c band on 
each sublattice is given by 



2T 
Jf 



(10) 



Here Q'^'^, Q^f and Q^'^ are the 2x2 sub-matrices of 
corresponding to the c, / and mixed sectors, respectively. 
The factors of two in the occupation numbers arise from 
summation over spin. 

In the slave-boson mean-field approximation, the orig- 
inal problem has thus been mapped onto a model of four 
bands with inequivalent hybridization \fJaboL between 
the Cq, and /„ bands. The resulting renormalized bands 
are defined by their dispersions i?km, = 1, 2, 3, 4, solu- 
tions of the eigenvalue problem 



Det^;-i(£:k,k) =0. 



(11) 



While it is not possible, in general, to obtain a closed 
form for the solution of (11), simplifications arise upon 



considering inter-sublattice (e.g. nearest-neighbor) hop- 
ping only, such that Ck = = 0. Further, if sublattice 
modulations of the hopping are absent, ek is real, and the 
c-electron kinetic energy is given by ±ek. In this case, 
one can invert the relation E = Emi^^s) and write 



J.bl - [X, + E)ip + E) 



X^ + E 

In particular, at the Fermi level, -E — 0, we have 



(12) 



ML = e(0) = 



A. 



(13) 



fiL may be interpreted as a c-band chemical potential ac- 
counting for the enlargement of the Fermi volume due to 



the Kondo effect, see Eq. (A3). In our case, 6i, 62 ^ at 
T = 0, and the resulting Fermi volume Vfl encompasses 
(2nc -I- 2) mod 2 electrons. 

In the mean-field equations, the restriction to inter- 
sublattice hopping only allows one to convert the mo- 
mentum summations into energy integrals 



dujpo (uj) . 



where ^0(1^) is the bare density of states (DOS) of the 
conduction electrons, 



po{io) = 1/AA^ [d{uj - ek) + S{uj + Ek) 



(14) 



The DOS is symmetric, po{i-o) = pq{—uj), due to Ck = 
^k — 0. Note that without sublattice modulations of 
the bare kinetic energy, po is identical to the bare local 
c-electron DOS on both sublattices A and B. 

In the following, the equations assume nearest- 
neighbor hopping only. To ease the numerical imple- 
mentation, the following calculations have been per- 
formed with either a semi-elliptic DOS, corresponding to 
a Bethe lattice with nearest-neighbor hopping, or a con- 
stant DOS, with no qualitative differences regarding the 
interplay between the two Kondo effects. In addition, we 
show some results for a square lattice, with and without 
next-neighbor hopping hopping, in Appendix [B} 

Finally, we note that non-local two-particle effects are 
not included in the present slave-boson approximation. 
In particular, there are no correlations between different 
local-moment spins, i.e., effects of RKKY interaction are 
neglected. At the mean-field level, this can be (partially) 
repaired by explicitly adding an inter-moment interac- 
tion which is then decoupled in the particle-particle or 
particle-hole channel. Such a procedure can, e.g., account 
for a spin-liquid phase competing with Kondo screening, 
but this is left for a future study. 



III. KONDO SCALES 

Let us start by recalling a few features of the slave- 
boson approximation applied to Kondo physics. The on- 
set of Kondo screening is in general signaled by an (ar- 
tificial) phase transition: For T > Tk, the mean field b 
vanishes while it is non-zero for T < Tk- For a single 
Kondo impurity coupled by an exchange coupling J to a 
bath with a DOS po{uj), the mean-field Kondo tempera- 
ture Tk{J) is the solution of 
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(15) 



For a constant DOS and in the weak-coupling limit, 
Tk{J) evaluates to 



Tk(J) ~ v/^2_^2exp [-l/(Jpo)] . (16) 
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Compared to the exact single-ion Kondo temperature, 
this mean-field Tk has the correct exponential de- 
pendence on J and po, and the prefactor matches the 
one-loop result for The artificial phase transition 

at Tk(J) will turn into the physically correct crossover 
upon including gauge fluctuations in the slave-particle 
description. 

Remarkably, the slave-boson approximation applied to 
a standard Kondo lattice, with bare c electron DOS po(i^) 
and Kondo coupling J, leads to an onset of screening 
at exactly the same Tk(J), simply because c electrons 
and local moments are decoupled for T > Tk(J). How- 
ever, coherent Fermi-liquid behavior may set in at a much 
lower temperature Tcoh (even at the mean- field level), due 
to protracted screening.l^Sl 

With two inequivalent ions, there will be two onset 
temperatures for screening, Ti and T2, i.e., upon increas- 
ing the temperature the mean field parameter b,^ will 
vanish at the temperature T^. From Ji > J2 we expect 
Ti > T2 . The regime T > Ti consists of conduction elec- 
tron decoupled from all local moments - this is a pertur- 
bative high- temperature regime. When T2 < T < Ti, the 
local moments on sublattice A are hybridized with the c 
electrons, while those on B are free ~ this is the regime of 
partial screening advertised in the introduction. Finally, 
below T2, all local moments are hybridized. As above, 
the coherence scale Tcoh, to be discussed in Sec. 
be significantly below T2. 



IV 
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A. Two-stage screening 

The Kondo scales Ti , T2 are readily obtained from the 
mean-field equations ^ in the limits 61 — > 0, 62 — J" 0, 
respectively. It is then easy to see that 



T^=Tk{Ji) 



(17) 



with T\^{J) in Eq. (15), because all moments are decou- 
pled for T > Ti. In contrast, T2 has to be determined 
from 



2 

T2 



(18) 



where pc2{^^T) is the local DOS of c electrons on sub- 
lattice B. While Pc2(w) is simply given by po(w -I- ^jl) for 
T > Ti, the onset of Kondo screening on the A sublat- 
tice renormalizes pc2 in a temperature-dependent fashion, 
which is captured by the following equation: 



Pc2(w,r) = a{uj,T)po(\uj + fj,\a{uj,T)). 



(19) 



Here, a{uj,T) is a "boost factor", whose general form is 
given in appendix [C] For T > r2 it reduces to 



a{uj,T) 



(20) 
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FIG. 3: Renormalized dispersions E(e) (thick line) and bare 
dispersions (thin line) as a function of e obtained from equa- 
tion ^ for (a) T2 < T <Ti and (b) T < T2. The bare c- 
electron DOS is constant, po = 1/(2D), with a half-bandwidth 
D = 10. 



Due to the non-trivial temperature and frequency de- 
pendence of Pc2, the screening scale T2 will not have the 
usual form ( 16 1 for the Kondo scale as a function of the 
Kondo coupling J2. Still, we expect that an increase of 
density of states at the Fermi level w = to enhance 
this scale. We have, omitting the explicit temperature 



dependence of the mean-field parameters 



Pc2(0) = 



pXi 



1 Po 




XiP 



(21) 



If /Oc2(0)/po(m) > 1? the density of states at the Fermi 
level is enhanced by the first Kondo screening and we 
expect, accordingly, a boost of the B-sublattice Kondo 
screening, i.e., T2[J2) > Tk{J2)- Alternatively, if 
Pc2{0) / Poip) < 1, the B-sublattice Kondo screening will 
be inhibited by the onset of the A-sublattice one. In gen- 
eral, the effect of a{u!, T) is two-fold: as a factor affecting 
the argument of po (uj) (which may induce a particle-hole 
asymmetry) and as an overall coefficient. For a constant 
bare density of states, used for most of the numerical 
calculations, the first effect is absent. Otherwise, it will 
induce quantitative changes depending on the shape of 
Pq - this happens, e.g., in the square-lattice case in Ap- 
pendix |b] 

Before discussing numerical results for the screening 
scales Ti^2, it is enlightening to look at the renormalized 
dispersions E{e) as obtained from Eq. ( |Tl| ) by solving the 
set of mean-field equations. Fig. [3] displays the generic 
resulting bands for T2 < T < Ti and T < T2 < Ti. For 
temperatures T2 <T < Ti (Fig. Isk) , the local moments 



6 




FIG. 4: Temperature evolution of the normalized density 
of states at the Fermi level pc2(0)/po(/i) (see Eq. (21 1) of c 
electrons in the second sublattice for J2 = at an electronic 
filling number (a) ric = 0.2, (b) ric — 0.9. The temperature 
scale is normalized by Ti, the first Kondo coupling is Ji = 8 
and the half-bandwidth is set to D = 10. The bare density 
of states po is taken constant in its support. Inset: averaged 
normalized density of states (Pc2(0))t/po(m) in ^ window of 
width ~ T around the Fermi level. 



of the B sublattice are unscreened; in the mean-field ap- 
proach, this is signaled by a flat band being located at 
the Fermi level w = 0. Hence, three renormalized bands 
form, one of which is narrow and primarily of / charac- 
ter. If only nearest-neighbor hopping is considered, the 
symmetry of the DOS imposes that at half-filling, = 1 
at = 0, this narrow band is completely flat and lo- 
cated aX u! — 0. Away from half-flUing, the width of 
this narrow band is then linear in \fi\. With an asym- 
metric DOS, the width of the narrow band still reaches 
its minimum at half-filling (but does not vanish). The 
dependence of this width on the chemical potential /i in- 
dicates a strong effect of the conduction-band filling Uc 
on the B-sublattice Kondo screening. Lowering the tem- 
perature further below T2 (Fig. [3]d) , a hybridization gap 
opens near the Fermi level, and the second /-level trans- 
forms into a weakly dispersive heavy band crossing the 
Fermi level. 



B. Competitive vs. co-operative Kondo screening 

We are now in the position to discuss whether the 
Kondo effects on the two sublattices compete or co- 



operate. To quantify this effect, we shall compare the 
screening-onset temperature T2 of the B-sublattice mo- 
ments in our model with Kondo couplings Ji > J2 to the 
screening-onset temperature of a homogeneous Kondo 
lattice with the same bare DOS where all Kondo cou- 
plings are set to J2 ~ the latter temperature is simply 
given by Tx{J2) according to Eq. (fTsl). 

As T2 is determined by Eq. (|18[), we first focus on 
the relevant B-sublattice local DOS pc2{<^,T). Fig. |4] 
shows the temperature evolution of the normalized DOS 
at the Fermi level Pc2(i^ = 0)/po(a*) for different values 
of the electronic filling ric ~ 0.2 and 0.9. Note that we 
have formally set J2 = 0, such that all temperatures 
obey T > T2- For both cases, Pc2(0) vanishes in an 
intermediate-temperature regime below Ti, as none of 
the quasiparticle bands crosses the Fermi level. The sys- 
tem is still metallic: the averaged density of states in a 
window of width ~ T around the Fermi level is finite, 
see the insets of Fig. [4j At very low temperatures, we 
recover a finite density of states /9c2(0), as we expect for 
a metallic state, except at quarter filling, Uc — 0.5, where 
Pc2 (0) = at all temperatures below Ti and a band insu- 
lator is obtained. This can be understood using a count- 
ing argument: at quarter-filling, and for J2 = 0, the large 
Fermi surface encompasses an even number -t- 1 = 2 
of particles and a band gap at the Fermi level is expected. 

The crucial difference between Figs.|4^ and|4)D is that 
for Uc = 0.2, and generically below quarter-filling, Pc2(0) 
is smaller than Po(m) ^-t low temperature. Accordingly, 
we expect the B-sublattice Kondo effect to be reduced 
as compared to the homogeneous case. In contrast, for 
ric — 0.9, and generically above quarter-filling, Pc2(0) is 
strongly enhanced. This enhancement reaches its maxi- 
mum at half-filling where the width of the narrow band 
reaches its minimal value as discussed earlier. Accord- 
ingly, we expect a sensitive boost of the second Kondo 
effect when J2 <C Ji . At quarter-filling, because the den- 
sity of states at the Fermi level vanishes at all temper- 
atures below Ti, a second Kondo screening occurs only 
above a finite critical value J2c- 

The B-sublattice screening scale T2, as obtained from 
Eq. (18 1, is displayed in Fig. [5] for different values of the 
conduction-band filling ric- At all fillings, the strength 
of the second Kondo scale is consistent with the expec- 
tations drawn from Fig. |4] for J2 ^ Ji : a reduction (en- 
hancement) of the B-sublattice c-electron local DOS re- 
sults in a suppression (promotion) of the B-sublattice 
Kondo effect. 

When ric ^ 1/2, the second Kondo effect is inhib- 
ited for all values of J2 < Ji (Fig. [5^), whereas it is 
promoted for 1/2 <C ric < 1 for all values of J2 < Ji 
(Fig. [sj;) . It is worth emphasizing that the suppression 
of T2 for small ric is a drastic effect: Even small dif- 
ferences between Ji and J2 can induce largely different 
Ti and T2, Fig. [5^. Close to quarter-filling, a crossover 
regime is found for which the B-sublattice Kondo effect 
is promoted for J2 ^ Ji and inhibited for intermediate 
values of J2 up to Ji (Fig. I5j3) . For these intermediate 
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FIG. 5: Evolution of the temperature scales r2(J2), rK(J2), 
where Tk{J2) is the Kondo scale obtained in the homogeneous 
Kondo lattice with a Kondo coupling J2, with the strength 
of the second coupling constant J2 < Ji for an electronic 
filling number: a) ric — 0.2, (b) — 0.6 and (c) ric = 0.9. 
The temperature scale is normalized by Ti , and the one of J2 
by Ji. The bare density of states po used is constant in its 
support. 



values of J2, it is not straightforward to deduce the effect 
of the A-sublattice Kondo effect on the B-sublattice one 
because the density of states at the Fermi level Pc2(0) 
is vanishing in an intermediate energy regime. One has 
then to consider the full integral in ( 18 ) instead of fo- 



cusing on its usually dominant u; — contribution. We 
found this transition regime to depend on the shape of 
the bare DOS. 



C. Kondo feedback 

So far, we have considered the effect of the A-sublatticc 
Kondo screening on the B-sublattice one, mainly in the 
extreme inhomogeneous case of small J2. We now discuss 




T/Ti 

FIG. 6: Temperature dependence of the A-sublattice con- 
densate 61 for J2 = 7, Ji = 8 at different filling numbers 
ric = 0.2,0.9. The temperature scale is normalized by Ti. 
The arrows point to the kinks occurring at the second Kondo 
scale T2. The bare density of states po used is constant in its 
support with a half-bandwidth D = 10. 



the "feedback" effect of the B-sublattice screening on the 
A-sublattice one. We have already emphasized that the 
A-sublattice Kondo scale is the same as the one obtained 
for the homogeneous case Ji — J2- That is because this 
scale depends on the bare density of states of c electrons. 
Fig. [6] shows the generic temperature evolution of the A- 
sublattice condensate bi at two different filling numbers 
Tic- A kink is observed as soon as the B-sublattice Kondo 
effect sets in. Below quarter-filling, when an inhibition 
of the B-sublattice Kondo screening is observed, bi is 
reduced below T2, whereas it increases above quarter- 
filling, saturating at low temperatures. This means that 
the two Kondo effects are either mutually co-operative 
(above quarter filling) or mutually competitive (below 
quarter filling) . Fig. [t] shows the generic temperature 
evolution of the average occupation number rid in the 




FIG. 7: Temperature evolution of the average occupation 
number rid in the first sublattice for J2 = 7, Ji = 8 at dif- 
ferent filling numbers ric ~ 0.2,0.9. The temperature scale 
is normalized by Ti and rid is normalized by the total filling 
number ric. The arrows point to the kinks occurring at the 
second Kondo scale T2. The bare density of states po used is 
constant in its support with a half-bandwidth D = 10. 
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A sublattice at two different fillings ric = 0.2,0.9. It 
turns out that for all fillings Uc this quantity increases 
as soon as the A-sublattice Kondo effect sets in, i.e there 
are more electrons on the A sublattice as compared to 
B (recall that, for T > Ti, both fillings are equal in our 
model). At the onset of the B-sublattice Kondo effect, 
Uci reduces again without necessarily recovering its high- 
temperature value. 

The fact that rid increases for all values of filling num- 
ber ric as soon as the A-sublattice Kondo effect sets in 
means that the latter sustains itself not only by involving 
low energy, i.e near the Fermi energy, conduction elec- 
trons in the second sublattice but also high energy ones. 
It shows also that the average number of c electrons in the 
B sublattice nc2 plays little role, compared to the den- 
sity of states Pc2 , in identifying a promotion or inhibition 
regime for the second Kondo effect. 



IV. LOW-TEMPERATURE PROPERTIES OF 
THE FERMI-LIQUID STATE 




T/T, 

FIG. 8: Total average entropy, normalized by 2 In 2,^^ at fill- 
ing number ric — 0.9 for different values of J2 — 0, 3, 7 and 
Ji — 8. The temperature scale is normalized by Ti. Tiny 
dashed lines are linear fit at very low temperatures. The ar- 
rows point to the kinks occurring at the second Kondo scale 
T2. The bare density of states po used is constant in its sup- 
port with a half-bandwidth D — 10. 



At sufficiently low temperatures, we expect our model 
system to behave as a Fermi liquid (neglecting possi- 
ble magnetic or superconducting instabilities and pro- 
vided that no Kondo insulator obtains). The Fermi liq- 
uid can be characterized by its linear-in-temperature spe- 
cific heat C{T)/T k, 7, constant magnetic susceptibility 
Xq and quadratic-in-temperature resistivity p{T) — pg w 
AT^ . This Fermi- liquid behavior will be realized below 
a so-called coherence scale Tcoh, which governs the low- 
temperature properties of thermodynamic and transport 
quantities and takes the role of an effective Fermi energy 
of the heavy quasiparticles. 

For the case of a homogeneous Kondo lattice, stud- 
ies using local self-energy approximations to the Kondo 
lattice model (i.e. dynamical mean-field theory and 
slave bosons) have shown that is different from Tk, but 
shares its exponential dependence on the Kondo coupling 
J JO-22 30 particular, their ratio has been found to be 
independent of J. In this section, we explore this issue 
in the case of two inequivalent local moments per unit 
cell. The effective four-band picture is particularly ap- 
propriate to describe the Fermi-liquid regime that sets in 
at very low temperatures. 

In the low-temperature limit, thermodynamic quanti- 
ties, such as the electronic entropy Se{T) or the electronic 
specific heat coefficient C iT) , are determined by the den- 
sity of states at the Fermi level p(0) 

2 2 

S,{T ^ 0) = ^p(0)T, C{T ^ 0) = '^p(0)T. (22) 
6 6 

Following Ref. [20l one may define a coherence scale 
in terms of the zero-temperature value of the average 
renormalized density of electronic states at the Fermi 
level p{Q) = \^^{pcv{^) + Pfvi'd)) as follows (recall 



^ 4p(o,r = o) ^^^^ 

so that the low-temperature expansion of thermodynamic 
and transport quantities may be written in terms of pow- 
ers of T/Tcoh- The same applies to our mean-field pa- 
rameters which are expected to saturate to their zero- 
temperature value sufficiently below Tcoh- As Tcoh is, 
strictly speaking, a scale characterizing the T = state 
only, we define Tcoh = rninlTcoh, Ti, T2}. (Note that 
Tcoh < T2 except in the proximity to Kondo-insulating 
states.) 



A. Entropy 

Fig. |8] shows the temperature dependence of the total 
average entropy for different coupling constants J2 above 
quarter-filling, here ric = 0.9. 

We first discuss the extremely inhomogeneous case 
J2 — 0, corresponding to decoupled B-sublattice mo- 
ments. Here, the entropy starts to drop slowly below Ti, 
remaining close to 4 In 2, then decreases steeply towards 
low temperature to reach 2 In 2 at T = 0. This residual 
entropy corresponds to the one of the unscreened local 
moments in the second sublattice.I^ Once this contribu- 
tion is subtracted, a linear dependence is found at low 
temperature consistent with the formation of a coherent 
Fermi liquid state, involving conduction electrons in both 
sublattices as well as the screened local moments on the 
A sublattice (while the B local moments are decoupled), 
with a finite "coherence" scale Tcoh. 



9 




FIG. 9: Evolution of the coherence scale Tcoh with the elec- 
tronic filling ric for different values of the coupling constant 
J2/J1 = 0,3/8,7/8,1. The bare density of states po used is 
constant in its support with a half-bandwidth D — 10. 



For finite values of J2, a second kink signals the onset 
of the second Kondo screening at T2 . Then a stronger de- 
crease is observed in the entropy below this scale. This 
decrease can be followed by a plateau-like feature close to 
half filling. (This feature is absent below quarter-filling.) 
At lower temperatures, the entropy is again linear in tem- 
perature signaling the emergence of a Fermi liquid regime 
out of the conduction electrons and the local moments in 
both sublattices. 



B. Coherence scale 

Fig. [9] shows the evolution of the coherence scale Tcoh 
with the filling number for different values of J2. In all 
cases, Tcoh approaches zero as ric —)■ as a result of 
"exhaustion" ^ 

Again, we start with the extremely inhomogeneous 
case J2 = 0. Here, a band insulator is obtained at 
quarter-filling and Tcoh is not defined. Above quarter- 
filling, the coherence scale for J2 — decreases down to 
zero at half-filling. This can be understood by looking 
at the renormalized band structure (see Fig. [s]). The 
width of the band crossing the Fermi level is of order /i, 
the chemical potential. As we increase the filling number 
Uc towards half-filling, this quantity, proportional to Tcoh, 
decreases towards zero. The argument for Tcoh going to 
zero at half-filling relies on the fact that the bare DOS 
is symmetric, and thus that fj, — at half-filling. How- 
ever, we observed that the decrease of Tcoh above quarter- 
filling holds even with the inclusion of next-nearest neigh- 
bor hopping, i.e. for a asymmetric bare DOS, albeit not 
towards zero at half-filling. 

For finite J2, the behavior of Tcoh is consistent with the 
observations regarding the promotion or inhibition of the 
B-sublattice Kondo effect, as its scale T2 sets an upper 
bound on Tcoh- This is more easily seen in Fig. 10 where 
we compare Tcoh to the homogeneous case's result. At 



T2(ii=l, Ji=8) 
fcoh(j2=7, Ji = 8) 
rcoh(j2=Jl=7) 



FIG. 10: Evolution of the coherence scale Tcoh and the sec- 
ond Kondo scale T2 with the electronic filling Uc for coupling 
constants J2/J1 = 7/8 (upper panel) and J2/J1 = 3/8 (lower 
panel). Also shown is Tcoh for the homogeneous case Ji = J2. 
The bare density of states pa used is constant in its support 
with a half-bandwidth D = 10. 



low fillings, the coherence scale Tcoh is smaller compared 
to the homogeneous case, consistent with the inhibition 
of the second Kondo effect: the mean-field parameters 
saturate at lower temperatures compared to the homo- 
geneous case. The opposite happens above, where Tcoh 
is increased, consistent with the promotion of the second 
Kondo effect. Its effect on Tcoh is most pronounced close 
to half filling. 

It is worth discussing the contributions to the total 
DOS p(0), which defines our Tcoh- Consider the ratio 
between the densities of states of the two sublattices: 



pM (i + if)pci(o) 



P2(0) 



M 

J 2b': 
A? 



(1 

1 (1 

a(0)2 (1 



-) Pc2(0) 



(24) 
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where a(0) is the boost factor given by 



q(0) = 



Jlbl 



1 - 



A2M 



see Eq. ( |C2[ ) in Appendix [C] for details. The ratio ([24]) 
is mainly controlled by two quantities: the boost factor 
a and the relative contributions coming from the heavy 
/-like bands forming below the Kondo scales. 

Fig. ( [TT| ) shows the contributions of each sublattice to 
the total DOS p{0) = l/(4f'coh) for different electronic 
fiUings He and couphngs J2. For < J2/J1 ^ 1, the 
total DOS is dominated mainly by contributions from 
the B sublattice for the lowest electronic filling up to 
Jic ^ 0.75 then decreases below the contribution from 
the A sublattice before converging to one half towards 
half filling. In this regime of Kondo couplings, and for 
low electronic filling, the main contribution comes from 
the /-like heavy band formed out of local moments in the 
B sublattice (regardless of the small boost factor). We 
note that the detailed behavior near ric ~ 0.8 depends on 
the input c-electron DOS and is different for the square- 
lattice case. 

For J2/J1 ^ 1, contributions from the A sublattice 
dominates for very low electronic fillings before being 
over passed by the contribution from the B sublattice 
already below quarter filling. The latter contribution at- 
tains its relative maximum at quarter filling before reduc- 
ing towards one half, both contributions being equal at 
half filling. Here the effect of the boost factor wins over 
the effect originating from the two /-like heavy bands 
and mainly sets the behavior of the ratio (24) 



An interesting feature for the different investigated 
shapes of DOS and < J2/J1 ratios, when only nearest- 
neighbor hopping is considered, is the convergence to- 
wards equal contributions from each sublattice when ap- 
proaching half filling, Uc — ^ 1. This is due to the sym- 
metry of the DOS that imposes /i = at half-filling. In 
this case, one can check from (24) that ^ — > causes 
pi{0) / P2{0) — > 1 for all values of J2 and Ji. 



V. DISCUSSION 
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FIG. 11: Contributions of each sublattice to the low- 
temperature average density of states at the Fermi level 
p(0) — l/4Tcoh, illustrated for different electronic filling Tie 
and J2/J1 = 3/8 (top) and J2/J1 = 7/8 (bottom), = 
Pcv + Pfv is the density of states of sublattice v. The bare 
density of states po used is constant in its support. 



effects are either mutually co-operative (ric > 1/2) or mu- 
tually competitive {ric < 1/2). A strong competition in 
the latter regime implies that small differences between 
Ji and J2 induce large differences between Ti and T2 as a 
result of non-trivial feedback effects, see Fig. [5^. There- 
fore, the regime of partial screening advocated in the in- 
troduction may be reached easily for small conduction 
band filling. 



A. Summary 



B. Relation to depleted Kondo lattices 



In this paper, we have investigated the interplay be- 
tween two sublattice Kondo effects in a bi-partite Kondo 
lattice model with two spin-1/2 local moments per unit 
cell, with Kondo couplings Ji > J2 on sublattices A,B. 
We have found that the conduction-band filling deter- 
mines how the B-sublattice Kondo effect is affected by the 
onset of A-sublattice Kondo screening: it is generically 
promoted for ric > 0.5 and inhibited for ric < 0.5. Our 
collected results thus show that quarter filling, Uc = 0.5, 
is the boundary separating regimes where the two Kondo 



Heavy fermion systems with randomly depleted local 
moments constitute an interesting route to interpolate 
between the physics of isolated Kondo impurities and 
coherent Kondo lattice behavior. Experimentally, deple- 
tion is realized by substituting magnetic by non-magnetic 
ions, as e.g. in Cci-ajLa^^Als. Theoretically, a number of 
papers have treated depleted Kondo lattices ,'25H23 ^j^h 
randomly distributed local moments with concentration 
rif < I, mainly using local self-energy approximations, 
combined either with CPA or with a numerically exact 
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treatment of spatial disorder. One of the remarkable find- 
ings is that nf = Tic marks a boundary between two dis- 
tinct regimes: For nf < ric the behavior is akin to that of 
dilute impurities, i.e., lattice effects appear to be small, 
whereas Uf > ric leads to Kondo lattice behavior. Al- 
though both regimes are Fermi-liquid-like at lowest tem- 
peratures, the crossover is rather pronounced e.g. in the 
coherence temperature Tcoh{nf) which is approximately 
constant for Uf < ric, but rapidly drops upon increasing 
Uf beyond Uc- The latter behavior may be related to "ex- 
haustion", i.e., the lack of conduction electrons to screen 
the moments. It has been speculated that the bound- 
ary = ric is inherited from the strong-coupling limit, 
J — >■ oo, where Uf — Uc represents a (random) Kondo 
insulator.'^ 

We believe that the two seemingly different problems 
of depleted Kondo lattices and Kondo lattices with in- 
equivalent local moments are related in the following way: 
For T2 < r < Ti in the latter problem, Kondo screening 
is "active" for half of the moments. This can be read 
as an n/ — 1/2 depleted Kondo lattice, albeit without 
randomness. Carrying over the results of Refs. [^5H?7l 
suggests that Uc < 1/2 then corresponds to a collective 
regime with exhaustion playing a role. As a result, con- 
duction electrons are strongly involved in screening A- 
sublattice moments, such that B-sublattice screening is 
suppressed and the two Kondo effects compete. In con- 
trast, for TT-c > 1/2 a single- ion picture is more appropri- 
ate: Kondo screening of A sites injects states near the 
Fermi level, such that B-sublattice screening is promoted 
and the two Kondo effects co-operate. These arguments, 
although somewhat simplistic, are consistent with our 
numerical results for pc2 in Sec. |III B[ 



uid with intra-cell dimer formation are possible, leading 
to a light metallic state with Fermi-liquid properties, sim- 
ilar to that in the bilayer Kondo lattice model of Ref . 1551 
However, exotic spin liquids may be formed e.g. from 
the B-sublattice moments only, such that a fractional- 
ized Fermi liquid may be realized in the regime of partial 
screening. A comprehensive study of these issues is left 
for future work. 



D. Experiments 

At present, relatively little detailed information is 
available about the low-temperature properties of the 
(RE)3Pd2oX6 compounds mentioned in the introduction. 
Clearly, a detailed characterization of the ordered phases, 
e.g., by neutron scattering, is desirable to put those in 
the global perspective advocated here. 

Based on our results, we propose to perform system- 
atic doping studies of these compounds, using dopants of 
different valence. As shown here, variations in the con- 
duction electron filling can be expected to significantly 
affect the different Kondo scales in the system. 
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Appendix A: Review of slave-boson results for the 
standard Kondo lattice 



C. Extensions 

A natural extension of the calculations presented here 
is to take into account non-local magnetic correlations, 
due to either RKKY interactions or a direct Heisenberg- 
like exchange. In addition to the standard Kondo lat- 
tice phenomenology, i.e., the occurrence of magnetic (or 
spin-liquid) phases for dominant inter-moment exchange, 
interesting physics can be expected in the regime of par- 
tial screening, see Fig. [l] For instance, for small J2 the 
B-sublattice moments are more susceptible to non-local 
effects than the A-sublattice moments, paving the way to 
either non-local B-sublattice "screening" or B-sublattice- 
dominated magnetic order. 

One might also speculate about variants of the Kondo- 
breakdown (KB) scenario for quantum criticality, orig- 
inally proposed as a possible description of some fas- 
cinating non-Fermi liquid properties in heavy fermion 
systems.l^^t2l] A central role in this scenario is played by 
so-called fractionalized Fermi-liquid states, formed out 
of conduction electrons weakly coupled to an exotic (i.e. 
fractionalized) spin liquid.'^ In the present situation with 
two local moments per unit cell, a "conventional" spin liq- 



Here, we recall some results from the mean-field ap- 
proximation of the Kondo lattice problem. An extensive 
presentation can be found in Ref. f3Ul, 

The mean-field equations are written as follows 



6 = vJ / du! pfc{(jj)nF{(jj), 
J —00 

/oo 
duj pf{LLi)nF{i^), (Al) 
-00 

/■oo 

Uc — 2 du! pc{io)nF{uj), 



where npd^) is the Fermi function and pc,Pf,Pfc are 
the spectral densities of the local single-particle Green's 
functions given as 



Pc(w) = 


TT 




Pfc{^^) = 


TT ■' 




Pf{uj) = 


TT ■' 


= (. + A)2^^(")' 
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where G~^{^) = G{i^ w + zO"'") is the retarded Green's 
function matrix and e(a;) = cj + /i — 
Notice that 



l)mod2 = 



(A3) 



2 
J 



^'^'^^ — ^ ^ tanh 



2Ik 



(A4) 



where ao = 1.13 is a numerical constant and 



duj 



Po(cj + m)-PoM 
2|w|po(m) 



(A5) 



(A6) 



T,, 



— X ' 

where A/i — fi^ — ^ and 



-i/Jpo 



i^o = exp 



Am 



(i'+p) 



pojuj + ^) - pajp) 



(A7) 



(A8) 



0.3 r 



0.2 



0.1 



can be understood as a chemical 



where /zl = /i — ^ 
potential corresponding to the enlarged Fermi volume, 
which encompasses n = ric + 1 electrons when b 0. 

The Kondo temperature Tk is defined as the temper- 
ature for which the effective hybridization oc b vanishes. 
It is obtained from the mean- field equations ( A2 1 in the 
limit 6, A 0. This gives 
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For /9o(m) 7^ Oj the main contribution to the integral in 
(A4) comes from a; = and any finite coupling J leads 



to Tk > 0. In the weak coupling regime J <^ D, one 
findPl 



In the Kondo lattice problem, the Kondo scale Tk is not 
necessarily the one that characterizes the behavior of the 
low-temperature regime. Generally, the thermodynamic 
quantities like the specific heat coefficient or the zero 
temperature magnetic susceptibility are characterized by 
a second energy scale given by 



From ( A5 1 and ( A7 ) , we see that the scales Tk and Tcoh 
have the same exponential dependence on the coupling 
constant J. They can differ considerably in their pref- 
actors. Their ratio depends generally on the electronic 
filling and the shape of the density of states. Well below 
^coh = minlTcoh, Tk}, the system behaves as a Fermi 
liquid for less than half-filling. 
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FIG. 12: Densities of states of the c-electrons in each sublat- 
tice, left for sublattice A and right for B, for T2 < T < Ti (top 
panels) and T < T2 (bottom panels) in the square lattice case 
with no next-nearest neighbor hopping. The bare density of 
states for T > Ti is shown in dashed line in each panel. The 
half bandwidth is D = 10. 



Restricting ourselves first to the case of nearest-neighbor 
hopping only, we can apply the numerical machinery as 
in the body of the paper, i.e., the momentum summa- 
tions can be converted into frequency integrals. To this 
end, one needs the c-electron density of states, which is 
given by 



TT^D 



-X 



1 



(Bl) 



where K,{x) is the elliptic integral of the first kind, and 
D — At where t is the nearest-neighbor hopping. This 
DOS is smooth and continuous, except at the step-like 




Appendix B: Results for the inhomogeneous square 
lattice 



In this appendix, we present some results for the 
Kondo lattice problem with two inequivalent local mo- 
ments, Eq. ([1]), on a square lattice of c-electron orbitals. 



FIG. 13: Evolution of the coherence scale fcoh and the Kondo 
scale T2 with the electronic filling ric for a coupling constant 
J2/J1 = 7/8 and comparison with the standard Kondo lattice 
result for the square lattice case with no next-nearest neighbor 
hopping. The half bandwidth is 73 = 10. 
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FIG. 14: Evolution of the temperature scales T2{J2) and 
Tk{J2), where Tk{J2) is the Kondo scale obtained in the 
homogeneous Kondo lattice with a Kondo coupling J2, with 
the strength of the second Kondo coupling J2 < Ji for an 
electronic filling — 0.2 and various values of the next- 
nearest hopping parameter t' . 



FIG. 15: Evolution of the temperature scales 72 (J2) and 
Tk{J2), where Tk{J2) is the Kondo scale obtained in the 
homogeneous Kondo lattice with a Kondo coupling J2, with 
the strength of the second Kondo coupling J2 < Ji for an 
electronic filling Uc — 0.9 and various values of the next- 
nearest hopping parameter t' . 



band edges and at the lo = van Hove singularity, where 
the DOS diverges as /Cn(w) ~ log '^^ 



Fig. 12 shows the temperature evolution of the local 
DOS for the c electrons in each sublattice for the square 
lattice with no next-nearest neighbor hopping. This evo- 
lution is consistent with Fig. |3] showing the renormalized 
dispersions at < T < Ti and T < Tj. 

Fig. [13] shows the dependence of the coherence and 
Kondo scales on the electronic filling number ric for 
J2 ~ Ji and J2 — Ji and how they compare to the coher- 
ence scale in the homogeneous case. This is qualitatively 
similar to Fig.jo] for low fillings, the coherence scale Tcoh 
is smaller compared to the homogeneous Kondo lattice 
case then bigger above quarter filling, when we approach 
half filling. Thus, the qualitative picture developed in the 
body of the paper is reproduced for the square lattice. 

Finally, we investigate the effect of the inclusion of 
next-nearest hopping, with a parameter t'. Non-zero 
t' results in non-vanishing momentum-dependent diag- 
onal terms = Ck = — 4i' cos (kx) cos (ky) in the kinetic 
matrix £k in Then, momentum summations in the 
mean-field equations (|9| can no longer be turned into an 
integration over energies. Instead, the mean- field equa- 
tions are now solved on a discrete lattice of 64^ sites. 

Figs. [M] and [15] show the dependence of the B- 
sublattice Kondo scale and the Kondo scale in the ho- 
mogeneous case on the electronic filling ric for various 
values of the next- nearest hopping parameter t' . 

Our conclusion on the promotion or the inhibition of 
the second Kondo effect is found to be robust against 
inclusion of next-nearest neighbor hopping, namely inhi- 
bition below quarter-filling and promotion above quarter- 
filling. 



Appendix C: Spectral densities 

The mean-field equations ^ can be written alterna- 
tively using the spectral densities of the local single- 
particle Green's functions 



(CI) 



K = V Jv I p/ciy(w)ni?(a;), 

— 00 

/•OO 

1 = 2/ duj pf^{uj)np{u!), 
2 / duj pcu{i^)np{uj), 



where np{uj) — 1/(1 -f e^l^') is the Fermi function. If we 
consider only nearest-neighbor hopping, we can have an- 
alytic expressions for the spectral densities pcv, Pff, Pfcu 



Pc^{uj) = a''{uj)po {e{uj)) , 



(C2) 



uj + Xi 



where 



' (A^ + uj)iJpbl - {Xp + uj) (/i + uj)) 
{Xf, + uj){J^bl - (A^ +uj){p + u)) 



defines the boost factor and t{uj) is given in Eq. ([12 
Here K^p are matrix elements of the 2x2 matrix 



A = 



1 

1 



(C3) 
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The equations (CI) and the relations between the dif- 



ferent spectral densities ( C2 ) for a given sublattice v are 



formally the same as in the homogeneous case Ji = J2 
(|A2l). 
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